Statics of non uniform Josephson junction parallel arrays: model vs. experiment 
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We study experimentally and numerically the zero-voltage supercurrent vs. magnetic field of 
non-uniform arrays of Josephson junctions parallel-connected by a superconducting stripline. The 
measured curves are complex, unique and in excellent agreement with numerical simulations using 
a specially developed model. Using this, we can optimize the arrays to have any desired interfer- 
ence pattern. Such new devices could find applications in magnetometry, quasiparticle mixers and 
detectors, flux-flow oscillators and superconducting electronics. 

PACS numbers: 74.81.Fa, 85.25.Dq 
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Josephson transmission lines have been experimentally 
and numerically investigated for their unique proper- 
ties deriving from the Josephson non-linearitjiSi^. Such 
systems include the "long" Josephson junction (LJJ) 
and the parallel array of "small" superconducting tun- 
nel junctions connected by a superconducting transmis- 
sion line, called discrete Josephson transmission line 
(DJTL]i^2.. "Small" vs. "long" qualifies the junction 
size relative to the Josephson length 2 



Xj = y/^ /{fi J c d) 
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where <f>o = h/2e is the reduced magnetic flux quantum, 
J c is the critical current density and d is the magnetic 
thickness^ between the junction electrodes. A DJTL 
can also be viewed as a parallel array of N — 1 SQUID 
\oo])^i±LLL. For magnetometry, such circuits have been 
considered recently because they are more sensitive than 
SQUIDs, the response of their total zero- voltage super- 
current /max to a magnetic field H contains much nar- 
rower bumps. I m ax{H) is then analog to the diffrac- 
tion pattern of N-slit optical gratings. Recently non- 
uniform junction distributions have been considered to 
synthesize an I ma x(-ff) curve with a single narrow peak 
at H = 0, allowing an absolute magnetic zero refer- 
ence as well as a higher sensitivity^. For DJTLs where 
the Josephson coupling is spatially discontinuous and ex- 
ists only in the junctions, the difference tp(x, t) between 
the superconducting phases of the top and bottom elec- 
trodes obeys an inhomogeneous sine-Gordon equation 20 . 
When the phases vary little from one junction to an- 
other, static and dynamic fluxons, i.e. 27r kinks of ip can 
occur in DJTLs. In practice, arrays of small junctions 
are more flexible than long junctions. They allow bet- 
ter impedance matching, their quasi-particle and surface 
dampings can be adjusted separately and as we will see 
here their interference patterns can be designed to suit 
specific needs. For microwave detection, small junctions 
have been successfully used to make quantum-noise lim- 
ited SIS heterodyne mixers and detectors from 100GHz 
to 1000GHz, primarily for radioastronomyiSiii. These 
devices use photo-assisted quasiparticle tunneling (the 
damping term in the sine-Gordon equation). Recently 



insulating layer (SiO) 
. 




base electrode (Nb) 



top electrode (Nb) 




FIG. 1: Photograph of circuit "A" (see Table 1). The lfim 2 
junctions are made out of a Nb/Al/AlOx/Nb trilayer, pat- 
terned using a self-aligned process, where the upper Nb layer 
is removed by RIE around the junction mesas and replaced 
by SiO. A Nb counter-electrode is then sputtered to contact 
junctions and make the Nb/ SiO/Nb rf tuning microstrip cir- 
cuit. This one, designed for broadband submillimeter-wave 
SIS mixing, is also a non- uniform DJTL or SQUIG. 



SIS mixer designs have mingled with DJTLs—^. In 
particular, it was demonstrated that a uniform DJTL 
circumvents the fundamental bandwidth limitation of 
resistive mixers^. We extended this concept to non- 
uniform arrays with various junction sizes and spacings, 
to optimize the impedance Z(f) over virtually any de- 
sired submillimeter-wave band (below the superconduc- 
tor's gap frequency) . Such junction arrays behave as mul- 
tipole bandpass filters 13 . They have also been recently 
proposed as submillimetric direct detectors In these ap- 
plications, the Josephson effect brings excess noise and 
instabilities. Therefore in practice, one quenches / max by 
applying in the plane of the oxide layer a magnetic field 
corresponding to $ = n<I>o, where usually n = 1 or 2. Yet 
the Josephson response to H of complex, non-uniform 
DJTLs is not straightforward, due to unknown flux in- 
ternal distributions, possible phase slippage, ac modes. 
This initiated our study of Josephson phenomena in SIS- 
array mixers. This Letter focuses on the static behavior. 
The non-uniform SlS-array mixers investigated (see 
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Table 1 and Fig. [TJ consist of N = 5 Nb/AlOx/Nb 
junctions embedded in a Nb/SiO/Nb stripline of width 
w = 5 — 7 ' firn^. The array length is I = 80 — 86/jm. 
Junctions have equal width and length (wj = lj w 1/xm) 
and their spacings were optimized so that .£(/) ~ 40Q 
over ~ 400 — 600GHz, for best rf coupling in that 
ranged. We successfully applied that design approach 
to SIS mixers for a spaceborne instrument 2 . The de- 
vices were made with a self- aligned, Selective Niobium 
Etch Process (SNEP) using optical lithography on fused 
quartz. Critical current densities J c range from 5 to 
l3kAcm~ 2 , junctions have V g = 2.85mV and subgap-to- 
normal resistance ratios of 15 — 20. The Swihart velocity 
was measured in LJJs located on the same wafer, to 
derive the specific junction capacitance for each batch. 
C s = 77 (resp. 81) ±3fF(xm- 2 for J c = 5 (resp. 10) 
kAcm~ 2 are close to the assumed value 80/F/ /im 2 . With 
a 250 — nm SiO layer, the Swihart velocity c in the ar- 
rays is ~ 0.14c. In the microwave problem, junctions are 
point-like and arrays are ID, since lj < w << A = c/f 
over the whole band. The "small junction" approxima- 
tion is justified because lj < Aj « 5/xm In a DJTL 
the size of fluxons at rest is defined by the discreteness 
parameter— Aj = Xj/a where a is the junction spac- 
ing. A ,/ > 1 implies that the flux permeates the ar- 
ray over several cells and that fluxons may move freely. 
For A, 7 < 1, fluxons must overcome an energy barrier 
to move from each cell to its neighbor. In our case, 
there is no periodicity a so our discreteness parameter 
is Aj = Aj/Z av where i av = l/(N — 1) is the average dis- 
tance between junctions. The design and performance of 
these non-uniform array mixers are detailed elsewhere 13 . 
Table 1. Non-uniform DJTL geometry 
(Dimensions of the junctions in fim, wj = lj = 1). 
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We model the device shown in Fig. [T] The phase differ- 
ence between the top and bottom superconducting layers 
satisfies in the static regime the following semilinear el- 
liptic partial differential equation^ 



g(x;y) 



A 2 ; 



sin ip = 0, 



(2) 



where g(x; y) is 1 in the Josephson junctions and out- 
side and we assume the same magnetic thickness d in all 
the device. This formulation guarantees the continuity 
of the normal gradient of if, the electrical current on the 
junction interface. 

The boundary conditions representing an external cur- 
rent input I and an magnetic field H applied along y are 

<Pv\l = [^Mo^] d /®o = Tvl/Ql) , 

= [H T Mo(l - d/$ a =HT(l-^, 

where < v < 1 gives the type of current feed. When 
v = 1 the current is only applied to the long boundaries 
y = 0, w (overlap feed) while v = corresponds to the 



inline feed. The quantities H and I are respectively the 
normalized magnetic field and current. 

We consider long and narrow strips containing a 
few small junctions of area lj x Wj placed on the line 
y = w/2 and centered on x — aj, j — l,n. For 
narrow strips w < 7r, only the first Fourier's mode 
needs to be taken into account 18 , then we write ip as 

<p(x,y) w vl{y - w/2) 2 /(lw) + cf>(x) . 
The first term of the equation takes care of the boundary 
conditions. Thus, we obtain the following equation for 4>: 

~<Pxx + G(x)[wj/(wX 2 j)\ sin0 = uj/l , 
where 7 = vl/w. We restrict ourselves to small junctions 
and we consider that the phase does not vary in the junc- 
tions. To keep a continuous model, we decrease the size 
of the junction without neglecting the current crossing 
it. Then q tends to a sum of delta functions 22,23 : 
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djS(x - a^ sin0 = v-, d 3 — 2 , 
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Despite its theoretical aspect, the delta function allows a 
detailed mathematical analysis of the solutions 2 ^. 

As opposed to standard models of this type of 
devic e) 1 ! 2 ' 6 , our approach does not neglect the variation 
of ip in the linear part of the circuit. This allows, for ex- 
ample, to show the differences between inline and overlap 
current feed. The main advantage of eq.® is that we can 
choose the position aj and size (wj,lj) of each junction 
and that we satisfy the matching conditions that exist in 
the original problem J2J . We have shown 2 ^ the periodic- 
ity of 7 m ax curve for an array where the junction positions 
are rational and analyzed the differences between inline 
and overlap current feed, the effect of one remote junc- 
tion on the array etc... There we established that for a 
device with sufficiently small junctions, the 7 max curve 
of eq.© tends to a simple function called the magnetic 
approximation. 

For the device A (resp. B) (see Table 1), dj ~ 
0.025. (resp. dj « 0.0178 .. .)• For such small coeffi- 
cients, the magnetic current which flows the junctions 
does not affect much the phase in the linear parts of 
eq.©. In this case, the phase gradient is proportional 
to the magnetic field H (see eq.@). We therefore as- 
sume that 4>(x) = Hx + c, so 7 = Y]j dj sm(Haj + 

c). To find the r y msx {H) curve for the magnetic ap- 
proximation, we take the derivative of 7 with respect 
to c. The values of c such that d-f/dc — are 



c m ax(H) = arctan 



Jmax (H) 



N 

3=1 



we obtain 



(4) 



The absolute value is to guarantee that this is a maxi- 
mum and not an extremum. This generalizes the stan- 
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FIG. 2: Measured and simulated (see eq.0) J m ax(#)in the 
V — state versus H at 4.2K for the arrays of junctions: A, B 
and F. Notice that for the device F, we simulate the device 
with one junction with an area ~ l.l^im 2 and the second 
» 0.9/im 2 . 



dard approach^. Fig|3 shows the measured and simu- 
lated supercurrent Imax(H) for the circuits A, B and C, 
with the H field applied along the y-axis. All circuits are 
from the same batch and F was used to calibrate the mea- 
surements and validate the numerical model. As an ad- 
ditional test of the model, we computed I max (H) curves 
of arrays with N = 5, 10 and 20 periodically spaced 
junctions, and obtained the same patterns reported in 
other works^. The experimental I m a, x (H) curves were 
all measured at A.2K . The magnetic field was produced 
by a superconducting electromagnet calibrated with a 
cryogenic Hall probe. Fig[5] is a blow-up of the region 
near H — where the features of the array geometry 
become visible. On a larger H scale (300 Gauss corre- 
sponds to <J> = <J>o in a lfim 2 junction) all curves fol- 
low the well-known Fraunhofer diffraction pattern given 
by the length lj of the junctions. Within this envelope, 
the Imax(-ff) curves of A and B follow complex yet re- 
producible patterns which are geometry-dependent, yet 
qualitatively similar. The / ma x(-ff) curves of A and B 
are undisputedly different, showing the extreme sensitiv- 




FIG. 3: Simulated Im&x(H) curve for an array of 5 junctions 
such (all lengths are given in jj,m) w = 5, a\ = 2, l\ = 21.2, 
h = 10.6, Is — 10.6, U = 21.2 and respective area of the 
junction are (given in /im 2 ): 0.24, 2.24, 4.96, 2.24, 0.24 and 
J c = 5kAcm~ 2 . 



ity of the Josephson "interference" to junction distribu- 
tion. In FigsEK an d B, the measured I max (H) curves 
were smaller than expected from I c . This can be due to 
the difference in magnetic thicknesses d in the junction 
and microstrip that we neglect^S. The simulated curves 
of A, B and F are homothetically adjusted on the current 
axis by respectively the factors 260, 270 and 201. For the 
H axis, the needed adjustment factors are 0.8, 0.8 and 
1 and in the case of F an offset of 2.8 Gauss is compen- 
sated for. For the SQUID (F) we had to slightly correct 
the area of one of the junctions to fit I max (H) and get its 
smooth behavior near the minima instead of the sharp 
behavior obtained for identical junctions^. The period- 
icity of imax(-ff) in the SQUID is H = 4.8 Gauss, close 
to the theoretical value Hq = 4.6 Gauss, derived from a 
calculated effective magnetic permeability ~ 1.68 in the 
Nb/ SiO/Nb stripline. The ratio of the bumps average 
widths in the N = 5 and N = 2 cases is ~ 8, which is 
the ratio of the circuits lengths. For / = 80/im, H w 0.6 
Gauss corresponds to $o- The narrow (0.5 — 0.8 Gauss) 
bumps for the N — 5 arrays therefore correspond to su- 
percurrent vortices extending over the whole circuit, in- 
dicating that the phase distributes itself over the whole 
device. 

Most experimental and numerical studies of fluxons in 
DJTLs have dealt with uniform arrays^iSiiL24 and large 
junction numbers (N > 10). In such a DJTLs with 
A j < 1, flux quanta experience pinning and therefore 
the "nuls" in I ma x(-ff) do not reach 1 = but a I m - m ^ 
floor which rises with Aj 1 ^. The non- uniform arrays re- 
ported here have different non-zero I m i n values for all 
minima and since all bumps in I max (H) are clearly seen, 
fluxon pinning is moderate. We expect mobile fluxons in 



the circuit A because A} w 1.78 > 1 (J c = bkAcm~ 2 ). 
We observed^ several constant- voltage branches in the 
IV Cs of A, B and other non-uniform arrays, due to 
Josephson resonances. Their frequencies /„ and their 
shapes suggest traveling fluxons. We have also seen such 
branches in the numerical solutions for given voltages. 

Fluxon propagation in LJJs has been exten- 
sively studied for submillimeter-wave flux-flow oscilla- 
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tors (FFOs) in superconducting integrated heterodyne 
receivers*^. Non-uniform DJTLs hosting fluxons offer 
advantages over LJJs for FFOs 14 . Whereas broadband 
rf coupling is a challenge for the very-low impedance 
LJJs, non-uniform DJTLs such as discussed in this 
Letter were optimally designed for this. Also, tunneling 
and stripline regions distinct in a DJTL allow to adjust 
separately quasiparticle loss and surface loss parameters 
which jointly determine the oscillator mode (resonant vs. 
flux- flow )£. Such non-uniform DJTLs could also be use- 
ful in fluxon-based electronics. 

Our mathematical model, validated by measures, al- 
lows to develop non-uniform DJTLs for magnetometrj^ 
(SQIFs). It can also be used to improve their perfor- 
mance as Josephson magnetic filters. For instance, mi- 
crowave circuits and sensors using quasiparticles should 
be most insensitive to H noise. To reduce Josephson 
noise in SIS mixers one biases them with a magnet on the 



zeroes of their Imax(U) curve. But mixer stability is lim- 
ited by how narrow the dips are at one zero. A specially 
designed non-uniform SIS array with zeroes considerably 
flattened out and set at specified H values would solve 
this problem. Just as an active antenna array can synthe- 
size a given microwave beam shape, non-uniform DJTLs 
can produce a given I milK (Il) curve and our model can be 
reversely applied to derive the required geometry. As an 
example, Fig[3| shows the simulated triangular-periodic 
^max(-ff) pattern corresponding to a N = 5 array with 
various junction spacings and sizes. This geometry can 
be realized by optical lithography. Other patterns where 
Im ax (H) is zero on a interval can also be generated. 

While at LERMA, M.-H. Chung of Taeduk Observa- 
tory developed the code for SIS mixing in DJTLs. We 
acknowledge Y. Delorme, F. Dauplay, B. Lecomte, A. 
Feret and J.-M. Krieg for their contribution to the SIS 
mixer work. 
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